Importance of interband transitions for the fractional quantum Hall effect in bilayer 

graphene 



Kyrylo Snizhko,^'^ Vadim Cheianov,^ and Steven H. Simon^ 

^Physics Department, Taras Shevchenko National University of Kyiv, Kyiv, 03022, Ukraine 

^Department of Physics, Lancaster University, Lancaster, LAI 4^B, UK 
^Rudolf Peierls Centre for Theoretical Physics, Oxford University, Oxford, 0X1 3NP, UK 

Several recent works have proposed that electron-electron interactions in bilayer graphene can 
be tuned with the help of external parameters, making it possible to stabilize different fractional 
quantum Hall states. In these prior works, phase diagrams were calculated based on a single Landau 
level approximation. We go beyond this approximation and investigate the influence of polarization 
effects and virtual interband transitions on the stability of fractional quantum Hall states in bilayer 
graphene. We find that for realistic values of the dielectric constant, the phase diagram is strongly 
modified by these effects. We illustrate this by evaluating the region of stability of the Pfaffian state. 



Recent experimental observation of the fractional steps 
in the Hall conductivity of ultraclean suspended graphene 
[H p and in graphene on the hexagonal boron nitride 
(h-BN) substrate |3] has opened a new chapter in the 
physics of the fractional quantum Hall effect (FQHE) 
where new states of matter might be observed. For ex- 
ample, the peculiarities of the single-particle spectrum of 
graphene are predicted to result in novel strongly corre- 
lated Hall fractions such as the SU(4) generalization of 
the Halperin-Laughlin state [4 . It has also been pro- 
posed that graphene allows for unprecedented tunabil- 
ity of the electron-electron interaction potential within a 
partially filled Landau level O [6] , allowing access to the 
experimental investigation of quantum phase transitions 
between different topological orders [7 and potentially 
stabilizing exotic topological states such as non-Abelian 
quantum Hall fluids [6^. Of particular interest from this 
perspective is bilayer graphene (BLG) in which the ma- 
trix elements of the Coulomb interaction within one Lan- 
dau level (also known as the Haldane pseudopotentials 
[8^ ) may be tuned in situ by the application of an exter- 
nal electric field breaking the symmetry between the two 
graphene layers [6^. 

The tunable quantum Hall effect in BLG was investi- 
gated theoretically in two recent papers |6| [7| where the 
stability conditions for several Abelian and non-Abelian 
FQHE states are found. Both groups use the same 
methodology based on the single Landau level approx- 
imation (SLLA). In this approximation the Hilbert space 
is constrained onto a given Landau level (LL) and the 
Hamiltonian takes the form of two-body interaction en- 
coded in Haldane pseudopotentials [8 . The SLLA Hamil- 
tonian is then exactly diagonalized for a small number 
of particles. The advantages of this approach are (a) 
its relative simplicity, and (b) its effectiveness in GaAs 
structures [9 due to the large values of cyclotron gaps 
achieved in the experiment and a large value of the di- 
electric constant, which suppresses transitions between 
different LLs [10 . However, as was discussed in [11 , 
BLG is a narrow gap semiconductor in the strong cou- 
pling regime, where the effects of virtual interband tran- 
sitions are essential. Moreover, as will be discussed in 



this work, the large cyclotron gap condition can be eas- 
ily violated due to the peculiar properties of the single 
particle spectrum. In this paper we develop a theoreti- 
cal framework incorporating such effects into the SLLA 
by means of perturbation theory. We are able to ana- 
lyze the conditions of the applicability of our approach, 
and we find that over a broad range of parameters these 
interband transitions can completely change the effec- 
tive interaction within this system. As an illustration of 
the importance of these effects, we use our approach to 
study the stability of the Moore-Read Pfaffian[12 state 
in the BLG system. We find dramatic differences from 
the phase diagram predicted in Ref. [6]. We conclude 
that all future analyses of the BLG system (or, indeed, 
all narrow gap semiconductors) will require more care- 
ful treatment of these effects than has been made in the 
past. 

The single-particle Hamiltonian for BLG in a perpen- 
dicular electric field is l3l 
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where ^ = ±1 for two different valleys, tt = + ipy^ 
Pi = —ihdi-\-eAi/c, A is a vector potential for the uniform 
perpendicular magnetic field 5, 2/7 is the single-particle 
mini-gap tunable by the electric field (for simplicity we 
consider U > 0), /i^ is the Bohr magneton and Sz is 
the z-projection of the electron spin. The Fermi velocity 
V ~ 10^ m/s, and the interlayer coupling is taken to be 
71 ^0.35 eV [H]. 

The eigenstates of the Hamiltonian ([T]) with |^| < 71 
are characterized by five quantum numbers: the valley 
index ^, the LL number n G Z+, the angular momentum 
jz = m + 1 with m G (Z+ — n), s = ±1 which denotes 
whether the state has positive or negative energy and 
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Sz = ±1/2. The wave functions in the n-th LL are: 



^n*'0n-l,m+l 
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where 'i\)nra is the nonrelativistic Landau wave func- 
tion 1^ in the n-th LL (negative LL wave functions 
are identicahy zero by definition). The amplitudes 
A^^ ,D^^ depend on U and B, and this depen- 

dence is crucial for tuning of the interaction matrix ele- 
ments. The single particle spectrum can be found from 
the equation 

((w-Cs)^ -7'n)((w+Ce)' -l\n-\)) = jHs'-u') (3) 

where, following [13 , we introduce Uc = 2v'^eB/'yic; we 
also introduce dimensionless parameters u = U/{hWc)^ 
7^ = ^i/ihujc), ^ = {E ^2fiBBSz)/{huJc), where E is the 
energy. 

The most quantitative method for the theoretical in- 
vestigation of the FQHE is exact diagonalization of the 
Hamiltonian of a small particle number system. However, 
this is possible only for relatively small Hilbert spaces, so 
going beyond the SLLA makes calculations practically 
impossible. The peculiarities of the BLG single-particle 
spectrum put rather tough constraints on the applica- 
bility of the SLLA. Figure [l^ shows the dependence of 
the several lowest LLs' energies on the magnetic field at 
U = 50 meV. Only the positive energy part of the spec- 
trum is shown. The negative energy LLs can be obtained 
by the particle-hole conjugation. We label each positive- 
energy LL by a pair of quantum numbers (n,^). One 
can see that at high magnetic fields the levels group into 
quasidegenerate doublets separated by energy of the or- 
der of hujc. Figure [TJd shows the dependence of the same 
LL energies on the gap parameter U at 5 = 10 T. Note 
that at large enough U (or small enough B) LLs cross. 
Thus the applicability of the SLLA puts an upper limit 
on U at any given B. At the same time, due to the 
small energy separation between the levels in each dou- 
blet (which is proportional to U) the SLLA is valid only 
for large enough U. 

Next we discuss the effect of the single-particle band 
structure on the interaction physics of FQHE in BLG. 
One can safely work in the approximation neglecting non- 
conservation of the spin and valley quantum numbers. 
In this approximation the following effects can be impor- 
tant. First, BLG is a narrow-gap semiconductor with 
strong vacuum polarization effects [11 . Second, the in- 
teraction can lead to the change of the order of filling of 
the LLs or to the appearance of the spin and valley unpo- 
larized states. Third, even when the dielectric constant 
is large virtual hopping between LLs can still renormalize 
the intra-LL interaction. 

Vacuum polarization. The virtual processes shown 
m fig. [2ji lead to a renormalization of the electron- 
electron interaction potential. The Fourier transform of 
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FIG. 1: (color online). The dependence of the lowest LLs' 
energies on (a) the magnetic field B at [/ = 50 meV, and (b) 
on the gap parameter [/ at B = 10 T. Each level is labeled 
by a pair of quantum numbers (n,^). 



the renormalized interaction potential can be represented 
as follows ^6 : 
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where v{q) = 2Tie^ /{qn) is the Fourier transform of the 
bare Coulomb potential, n is the dielectric constant of 
the environment [17 and n((7,cj) is the polarization func- 
tion. The approximation neglecting the retardation ef- 
fects {uo = 0) is justified as long as we are interested 
in energies much smaller than the inter-LL gaps. We 
take the polarization function from the random-phase 
approximation calculation for BLG with magnetic field, 
which can be justified within the 1/N expansion [11] 
(TV = 2 spin projections x 2 valleys = 4). Since Il{q^uj = 
O) (X. q^ screening is not efficient at large distances, how- 
ever it strongly affects the first few Haldane pseudopo- 
tentials (corresponding to distances of the order of mag- 
netic length) which have the most significant impact on 
the stability of any FQHE state. 

Population reversal of LLs. The order of levels in Fig. 
[ijprescribes the natural order of filling of the LLs by elec- 
trons in the independent electrons approximation. How- 
ever, we find that for integer filling fractions the electron- 
electron interaction leads to a reversal of this natural or- 
der in a significant part of parameter space. For example, 
the Coulomb energy of the fully filled (2, +1) LL is less 
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(a) (b) 

FIG. 2: Feynman diagrams [18] showing renormalization of 
the electron-electron interaction due to (a) the vacuum po- 
larization processes, and (b) the simplest processes involving 
virtual hopping of one or both of the two interacting electrons 
from the n-th LL to the n^-th LL. 



than the one of the fully filled (2,-1) for U > 0. In 
the region where the interaction is strong compared to 
the gap between the two levels this leads the fully filled 
(2, +1) LL having lower total energy than the fully filled 
(2, —1) level. Thus the former will be filled before the 
latter. When the quasidegenerate levels are from dif- 
ferent valleys, valley-unpolarized state can be preferred, 
particularly for fractional filling. Furthermore, when the 
quasidegenerate levels are from the same valley (as in 
the n = and n = 1 case) level mixing can occur. 
These are interesting effects which are, however, beyond 
the scope of this paper. Whereas for fractional filling 
such effects are much more difficult to analyze, popula- 
tion reversal at the integer filling fraction is an indicator 
of strong violation of the SLLA. Thus we constrain our 
analysis to the region of the parameter space where no 
population reversal occurs at integer filling. This guar- 
antees that the states are valley polarized. For the valley 
polarized states one can investigate whether the state 
is spin-polarized. Generally, spin-unpolarized states are 
not favored by Coulomb repulsion unless the potential is 
hollow-core, which we find to be not true in our case. 

Renormalization of pseudopotentials due to virtual 
hopping. The SLLA is exact in the limit of infinite energy 
difference AE between LLs. For finite AE^ the SLLA 
pseudopotentials acquire corrections due to virtual tran- 
sitions between the LLs such as, for example, shown in 
Fig. [2]3[19 . Such corrections are theoretically tractable 
only in the perturbative regime (when they are small), 
however even the presence of small corrections may dra- 
matically affect the phase diagram due to the extreme 
sensitivity of the FQHE states to the details of the inter- 
action. In this work we take the virtual hopping correc- 
tions into account in second order perturbation theory 
(Fig. [2JD). We restrict the region of validity of our con- 
sideration by requiring the third order corrections to be 
smaller than the second order ones. 

To investigate the role of these effects on the stability 
of FQHE states we focus on the Pfaffian state. Our choice 
is motivated by the following considerations. First, this 
state is particularly sensitive to the details of the inter- 
action so it is a good illustration for our analysis. Sec- 
ond, the stability of this state in BLG was investigated 
in Refs. [5, 6 in the SLLA approximation but without 
these effects taken into account, so we can compare the 
phase diagrams. Third, the Pfaffian itself is an important 
state because it is an example of the non-Abelian topo- 



logical fluid. The tunable parameters are the magnetic 
fleld the electric field which determines the mini-gap 
parameter U and the dielectric constant which con- 
trols the deviation from the naive SLLA (which is ex- 
act for K: 00). We can also choose the half-filled LL 
number. Here we will concentrate only on the two lev- 
els: (1,-1) and (2,-1). The (1,-1) level wave func- 
tion is constructed from the nonrelativistic n = and 
n = 1 LL wave functions, the (2,-1) level wave func- 
tion is constructed from the nonrelativistic n = 0, 1, 2 LL 
wave functions. In both cases one can tune the pseudopo- 
tentials close to their values at the nonrelativistic n = 1 
LL, where the 5/2 state is observed in GaAs. 

The tuning mechanisms are, however, different for the 
two levels. Amplitudes of the wave function (|2| in the 
(1,-1) LL show little dependence on U so the main con- 
trol parameter is B. In contrast, the amplitudes of the 
wave function in the (2,-1) LL mainly depend on one 
parameter which is the U/hWc ratio, so both B and U 
can be used for tuning. 

The main factors determining deviation from the naive 
SLLA for the two levels are the polarization and virtual 
hopping to the nearby levels. For the (1,-1) LL this is 
hopping to the (0, —1) and the (2, —1) LLs, while for the 
(2,-1) LL the important hopping is to the (3,-1) LL. 
In addition to this, for the (2, —1) LL, it is important to 
consider effects of mixing with the (0,-1) and possible 
population order reversal with (2,+l). The latter are 
important factors restricting the region of applicability of 
perturbative analysis, however, when suppressed they do 
not lead to a renormalization of the intra LL interaction. 

Figures [3^ and|3]3 show the regions of the applicability 
of perturbative analysis for different values of n for the 
(1,-1) and the (2,-1) LLs respectively. For the (1,-1) 
LL the region is bounded from above by the condition 
of small hopping to the (2,-1) LL, the lower bound is 
due to the condition of small hopping to the (0, —1) LL. 
At small enough magnetic fields, at least one of the con- 
ditions is violated at all values of U . For the (2, —1) 
LL the region's upper bound is due to the condition of 
small mixing with the (0,-1) LL, while the lower and 
the left bounds are due to the condition of absence of the 
population reversal. The thick black line shows where the 
maximum overlap with the Pfaffian for the bare Goulomb 
interaction is achieved. One can see that for small dielec- 
tric constants this line lies outside the region of validity 
of perturbative analysis; however, for large enough n they 
intersect near = 50 meV in both cases. This happens 
at > 10 and A. > 6 for the (1, -1) LL and the (2, -1) 
LL respectively. 

Figures and show the dependence of the overlap 
of the exact ground state of the system with the Pfaf- 
fian on the magnetic field and the dielectric constant at 
/7 = 50 meV for the (1,-1) and the (2,-1) LLs respec- 
tively. The region where perturbative analysis is not ap- 
plicable is hatched. As one can see, for the (1,-1) level, 
a high overlap up to 0.94 (compare with non-relativist ic 
n = 1 level overlap of 0.7) is achieved for all admissible 
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FIG. 3: (color online). The region ((a) - for the (1,-1) LL, 
(b) - for the (2,-1) LL) of the applicability of perturbative 
analysis for fixed values of k = 5, 10 and 15. The size of the 
region increases with increasing k. The thick black line shows 
where the maximum overlap with the Moore-Read Pfaffian 
for the bare Coulomb interaction is achieved. 
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FIG. 4: (color online). Color plot of the overlap of the ground 
state with the Moore- Read Pfaffian for 12 particles at [/ = 
50 meV as a function of the magnetic field B and the dielectric 
constant k ((a) - for the (1, -1) LL, (b) - for the (2, -1) LL). 
Contours show the lines of constant overlap. The region where 
perturbative analysis is not applicable is hatched. 

values of Hi near B = 8 T. For the (2,-1) level a high 
overlap up to 0.93 is achieved near 5 = 20 T, also for 
all admissible values of k. However, the behavior of the 
high-overlap region is different for the two cases. The re- 
gion is situated at the same magnetic field but becomes 
narrower with decreasing for the (1,-1) LL. For the 
(2,-1) LL the region also gets narrower with decreasing 
hz but its position shifts to higher values of the magnetic 
field. 

The authors of [6 found that in the (1,-1) LL, high 
overlap is achieved in the region near B = 10 T. We 



find that the region of high overlap is situated there for 
large enough values of hz. However, for smaller values 
of Ai: < 10 the effect of level mixing becomes significant 
which makes observation of the Pfaffian state unlikely. 
The (2,-1) LL was also considered in [6 where it was 
concluded that the maximal overlap with the Pfaffian on 
this level is less than 0.6. Our results do not support this 
conclusion (even for the bare Coulomb interaction). 

The previous consideration shows that BLG can be 
tuned into the regime of high overlap with the Pfaffian; 
however, one needs a higher dielectric constant than the 
usual hz ^ 2.5 for graphene on a Si02 substrate. This is 
experimentally achievable. For example, on a Hf02 sub- 
strate [20l Hi is around 12.5. For this value of = 12.5 the 
required magnetic fields (around B = lOT for the (1, —1) 
LL and around B = 20T for the (2,-1) LL) are quite re- 
alistic. The gap to the first excited state at these param- 
eter values is around 2 K and 8 K for the (1,-1) and the 
(2, —1) LLs, respectively. With increasing magnetic field 
we find that the gap monotonically increases to the val- 
ues of around 20 K and 17 K, respectively, at 5 = 30 T. 
At the same time the overlap decreases to around 0.5, 
which is still fairly large. [21] This result, obtained for a 
finite number of particles, suggests that the system may 
still be in the same topological phase at higher magnetic 
fields. 

To conclude, we have analyzed the influence of inter- 
Landau level transitions on the phase diagram of the 
FQHE states. We find that the SLLA can only be used 
under quite stringent restrictions. For moderate values of 
the dielectric constant these effects can be taken into ac- 
count perturbatively leading to some modification of the 
phase diagram, with dramatic modification outside the 
perturbative approach applicability region where Pfaffian 
observation is unlikely. 

We thank Prof. V. Fal'ko for useful discussions. We 
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onalization, especially N. Regnault. We also thank the 
authors of the JaxoDraw program, which was used to 
draw the Feynman diagrams. 
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In this section we provide additional details regarding 
the stability of the Pfaffian state. 

Figures la and lb show the gap between the exact 
ground state of the system with 12 particles and its ex- 
act first excited state. In these figures, the gap is plotted 
as a function of the magnetic field and the dielectric con- 
stant at /7 = 50 meV for the (1,-1) LL (Fig. la) and 
the (2, —1) LL (Fig. lb). The region where perturbative 
analysis is not applicable is hatched. Figures 2a and 2b 
show the same gap in the units of typical Co ulomb e nergy 

/ where I is the magnetic length I = y^hc/eB. 

As one can see, for both levels gaps increase as the mag- 
netic field increases (staring from points where maximum 
overlap occurs) and as the dielectric constant decreases. 
This can be partially (for magnetic field) or fully (for the 
dielectric constant) attributed to the increase of /ni. 
This result, obtained for a finite number of particles, sug- 
gests that the system may still be in the same topological 
phase at higher magnetic fields. The same holds for low 
dielectric constants down to boundary of the perturbative 
approach applicability region: for smaller dielectric con- 
stants the physics is essentially non-single Landau level. 

Figures 3a and 3b show the dependence of the gap on 
the magnetic field at = 50 meV and n = 12.5 for the 
(1,-1) and the (2,-1) LLs respectively. Figures 4a and 
4b show the same data plotted in the units of j ni. One 
peculiarity may be worth noticing: the gap has a local 
minimum near the maximum overlap point. 

Figures 5a and 5b show the dependence of the gap on 
the dielectric constant at /7 = 50 meV and fixed B for 
the (1,-1) and the (2,-1) LLs respectively (5 = 8 T 
for (1,-1) LL and 5 = 19 T for (2,-1) LL). Figures 
6a and 6b show the same data plotted in the units of 

I kI. This illustrates once more that the gaps increase 
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FIG. 1: Color plot of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV as a function of the magnetic field B and the dielec- 
tric constant n ((a) - for the (1, —1) LL, (6) - for the (2, —1) 
LL) . The region where perturbative analysis is not applicable 
is hatched. 
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FIG. 2: Color plot of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV as a function of the magnetic field B and the dielec- 
tric constant k ((a) - for the (1, —1) LL, (6) - for the (2, — 1) 
LL) . The region where perturbative analysis is not applicable 
is hatched. 
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FIG. 3: Dependence of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV and k — 12.5 as a function of the magnetic field B 
((a) - for the (1, -1) LL, ih) - for the (2, -1) LL). Only the 
part where perturbative analysis is applicable is shown. 



as the magnetic field increases and as the dielectric con- 
stant decreases, however, the "intrinsic" gaps (gaps in 
units of / nt) increase as the magnetic field increases 
but decrease as the dielectric constant decreases. This 
shows the tendency towards the loss of stability by the 
Pfaffian state at small dielectric constants. 

Figure 7 compares the data on the overlap with Pfaf- 
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FIG. 4: Dependence of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV and hi — 12.5 as a function of the magnetic field B 
((a) - for the (1, -1) LL, (b) - for the (2, -1) LL). Only the 
part where perturbative analysis is applicable is shown. 
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FIG. 5: Dependence of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV and fixed B as a function of the dielectric constant k 
((a) - for the (1, -1) LL at B = 8 T, (b) - for the (2, -1) LL 
at B = 19 T). Only the part where perturbative analysis is 
applicable is shown. 



FIG. 8: Dependence of the overlap with Pfaffian and gap to 
the first excited state on the number of particles N (N — 
8,10,12,14) for (1,-1) LL at [/ = 50 meV, n = 12.5 and 
B = 8T (black), B = 20 T (blue), 5 = 30 T (purple). The 
gap is presented in units of typical Coulomb energy e'^/ni. 




FIG. 6: Dependence of the gap between the ground state 
and the first excited state computed for 12 particles at [/ = 
50 meV and fixed B as a function of the dielectric constant k 
((a) - for the (1, -1) LL at B = 8 T, (b) - for the (2, -1) LL 
at S = 19 T). Only the part where perturbative analysis is 
applicable is shown. 
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fian and the gap to the first excited state for the (1,-1) 
and the (2,-1) LLs in BLG and the non-relativistic n = 1 
LL. The data are shown for different numbers of particles 
N = 8, 10, 12, 14. For the (1, -1) and the (2, -1) LLs in 
BLG external parameters are set to be near the maxi- 
mum overlap region. The gap is presented in units of the 
typical Coulomb energy e'^/hzi. The results for BLG lev- 
els appear to be more stable as the number of particles 
grows than for the non-relativistic system. 

Figures 8, 9 present the dependence of the overlap and 
the gap at several parameter points for (1,-1) LL and 
(2,-1) LL respectively. Points are taken to be at the 
same mini-gap and the dielectric constant values as in 
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FIG. 7: Dependence of the overlap with Pfaffian and gap to 
the first excited state on the number of particles N (N = 
8, 10, 12, 14). Solid line is for non-relativistic n = 1 LL. 
Dashed line is for (1,-1) LL at U — 50 meV, k = 12.5 and 
B = 8T. Dot-dashed line is for (2, -1) LL at [/ = 50 meV, 
K = 12.5 and B = 19 T. The gap is presented in units of 
typical Coulomb energy e'^/hvi. 



FIG. 9: Dependence of the overlap with Pfaffian and gap to 
the first excited state on the number of particles N (N — 
8,10,12,14) for (2,-1) LL at = 50 meV, k = 12.5 and 
B = 19 T (black), B = 13 T (dark green), 5 = 15 T (green), 
B = 30 T (purple). The gap is presented in units of typical 
Coulomb energy / 



the previous figure, but the magnetic field changes. For 
(1,-1) LL we took the magnetic fields to be only higher 
than the maximum overlap region (since lower fields are 
outside the region of applicability of perturbative ap- 
proach). For (2,-1) LL we took both higher and lower 
magnetic fields. As one can see the overlap at the max- 
imum point seems to be more stable as system size in- 
creases than at the other points. However, the overlap at 
higher magnetic field for (2,-1) LL is still quite high and 
quite stable. The situation with the gaps appears sim- 
ilar. Though this does not give any conclusive answer, 
one might suggest that for (1,-1) LL the Pfaffian state is 
stable in the neighbourhood of the maximum overlap re- 
gion, while for (2,-1) LL it may be stable also at higher 
magnetic fields. 



